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We study B — > K*^*^ t^l~ decays based on a large data sample of 657 million BB pairs collected 
with the Belle detector at the KEKB e"'"e~ collider. The differential branching fraction, the isospin 
asymmetry, the K* polarization, and the forward-backward asymmetry {Afb) as functions of 
are reported. The fitted Afb spectrum tends to be shifted toward the positive side from the 
SM expectation. The measured branching fractions and lepton flavor ratios (electron/muon) are 
B{B K*£+r) = (10.8li;J ± 0.9) x 10^^ B{B Kl+T) = {A.9,t°ol ± 0.3) x 10"^ Rk* = 
1.21±0.25 ± 0.08, and Rk = 0.97±0.18 ± 0.06, respectively. 

PACS numbers: 13.25 Hw, 13.20 He 



The b — > st^ transition is a flavor-changing neutral 
current (FCNC) process, which, in the Standard Model 
(SM), proceeds via either a Z/7 penguin or a box di- 
agram at lowest order. The efi^ective Wilson coefficients 
C7, Cg, and Cio describe the amplitudes from the electro- 



magnetic penguin, the vector electroweak, and the axial- 
vector electroweak contributions, respectively. These am- 
plitudes may interfere with the contributions from non- 
SM particles Therefore, the branching fraction and 
the lepton forward-backward asymmetry (^Aps) as func- 
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tions of dilcpton invariant mass in 6 — s- s£~^£~ provide 
information on the coefficients associated with certain 
theoretical models 0, and are also sensitive probes for 
the presence of new physics. 

In this paper, we present measurements of the dif- 
ferential branching fractions and the isospin asymme- 
tries as functions of = M'^^c^ for B K*£'^£~ and 
B K£~^£~ decays. The K* polarization and Afb for 
B K*£'^£~ decays as functions of are presented as 
well. A data sample of 657 million BB pairs collected 
with the Belle detector at the KEKB e+e~ collider is ex- 
amined. Charge-conjugate decays are implied through- 
out the paper. Equal production of B'^B^ and B^B^ 
pairs are assumed throughout this paper. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer located at the KEKB collider 0], and consists 
of a silicon vertex detector (SVD), a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid that provides 
a 1.5 T magnetic field. An iron flux-return located out- 
side the coil is instrumented to detect mesons and 
to identify muons (KLM). The detector is described in 
detail elsewhere 

We reconstruct B K^*H'^£~ signal events in 10 final 
states: K+n-, KItt+ , K+tt°, K+ , and Kg for K^*\ and 
combine with either electron or muon pairs. All charged 
tracks other than the Kg tt+tt^ daughters are re- 
quired to have a maximum distance to the interaction 
point (IP) of 5 cm along the beam direction (z) and 0.5 
cm in the transverse plane (r-</)) . A track is identified as 
a K'^ (tt^) if the kaon likelihood ratio is greater (less) 
than 0.6 (0.4); the kaon likelihood ratio is defined by 
TZk = Ck/{^k + Ctt), where Ck {^n) denotes a like- 
lihood that combines measurements from the ACC, the 
TOF, and dE/dx from the CDC for the K+ (7r+) tracks. 
This selection is about 85% (89%) efficient for kaons (pi- 
ous) while removing about 97% (91%) of pious (kaons). 
In addition to the information included in the kaon like- 
lihood ratio, muon (electron) candidates are required to 
be associated with KLM detector hits (ECL calorimeter 
showers). We define the likelihood ratio TZ^ {x denotes /i 
or e) as TZ^ = Cx/{Cx + Cnot-x), where and Cnot-x 
are the likelihood measurements from the relevant detec- 
tors 01. We select ^± candidates with TZf, > 0.9 (0.97) 
if > 1 GeV/c (0.7< <1.0 GeV/c). These require- 
ments retain about 80% of muons while removing 98.5% 
of pions. Electron candidates are required to have TZ^ > 
0.9, TZfj, < 0.8, and pe > 0.4 GeV/c. These requirements 
retain about 90% of electrons while removing 99.7% of 
pions. Bremsstrahlung photons emitted by the electrons 
are recovered by adding neutral clusters found within a 
50 mrad cone along the electron direction. The energy 
of the additional photon is required to be less than 0.5 



GeV. 

Two oppositely-charged tracks are used to reconstruct 
Kg 7r+7r~ candidates. The invariant mass is required 
to be within the range 483-513 MeV/c^ (±5 times the K^ 
reconstructed-mass resolution). Other selection criteria 
are mainly based on the distance and the direction of the 
Kg vertex and the distance of daughter tracks to the IP. 
For tt" 77 candidates, a minimum photon energy of 50 
MeV in lab frame is required and the invariant mass must 
be in the range 115 < M^-y < 152 MeV/c^ (±3 times the 
TT*^ reconstructed-mass resolution). Requirements on the 
photon energy asymmetry, \E}^ — E^\/{E}^ + E^) < 0.9, 
and the minimum momentum of the tt" candidate in the 
lab frame, Pj^o > 200 McV/c, suppress the combinatorial 
background. 

i?-meson candidates are reconstructed by combining 
a K^*^ candidate and a pair of oppositely charged lep- 
tons, and selected using the beam-energy constrained 
mass A/bc = \/ ^beam ~ Pb energy difference 

AE = Eb — Eheum where Eb and pB are the recon- 
structed energy and momentum of the B candidate in 
the T(4S') rest frame and ii'beam is the beam energy 
in this frame. Bremsstrahlung photons are included in 
the calculation of the momenta of electrons and hence 
are included in the calculations of Mbc, -Ebeam and q^. 
We require B-meson candidates to be within the region 
Mbc > 5.20 GeV/c2 and -35 (-55) < A£; < 35 MeV 
for the muon (electron) modes. The signal region is de- 
fined by 5.27 < Mbc < 5.29 GeV/c^. For K* modes, 
the Mk-k candidate (signal) region is defined as Mktt < 
1.2 GeV/c2 {\Mk^ - rriK' \ < 80 MeV/c^). 

The main backgrounds are continuum e+e^ qq 
{q — ti, d, c, s) and semileptonic B events. A Fisher 
discriminant including 16 modified Fox-Wolfram mo- 
ments is used to exploit the differences between the 
event shapes for continuum qq production (jet-like) and 
for BB decay (spherical) in the e+e~ rest frame. We 
combine 1) the Fisher discriminant. 2) the missing mass 
A/niiss = V^^miss "-Pmiss' 3) the angle between the mo- 
mentum vector of the reconstructed B candidate and the 
beam direction (cos 0b), and 4) the distance in the z 
direction between the candidate B vertex and a vertex 
position formed by the charged tracks that are not as- 
sociated with the candidate B-meson into a single likeli- 
hood ratio TZ = Cs/{Cs + Cqq), where Cs {Cqq) denotes 
the signal (continuum) likelihood. For the suppression 
of semileptonic B decays, we combine the Fisher dis- 
criminant, Mmissj cos 9b, and the lepton separation near 
the IP in the z direction to form the likelihood ratio 
TZb = Cs/{Cs + ^bb)^ where C^-g is the likcfihood for 
semileptonic B decays. 

Combinatorial background suppression is improved by 
including q^ and _B-flavor tagging information [7| , which 
is parameterized by a discrete variable gtag indicating the 
fiavor of the tagging i?-meson candidate and a quality 
parameter r (ranging from for no flavor information to 
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1 for unambiguous flavor assignment). Selection criteria 
for TZ and TZb are determined by maximizing the value of 
S/ y/S + B, where S and B denote the expected yields of 
signal and background events in the signal region, respec- 
tively, in different (g^, grec • 9tag • r) regions, where 
the charge of the reconstructed B candidate. Events with 
Qicc ■ Qtag • ''' close to —1 are considered to be well tagged 
and are unlikely to be from continuum processes. For the 
Kg£^£~ modes, only the dependence on r is considered. 

The dominant peaking backgrounds are from the B 
J/tpX and ip'X decays and are rejected in the regions 
(the limits are given in units of GeV^/c^): 

8.68 < q^in+fi-} < 10.09 , 

12.86 < q^{^l+^l-) < 14.18 , 

8.11 < g2(e+e-) < 10.03 , 

12.15 < q^{e+e-) < 14.11 . 

The decay B+ J/tlj{i'')h+ {h+ = K+,tt+) can also 
contribute to the 5+ K^tt^ i^i^ ij.^ and Kgir^n'^fi^ 
samples if a muon from J/'ip[tp') is misidentified as a 
pion and another non-muon track is at the same time 
misidentified as a muon. We remove such events from 
the two samples with the requirement —0.10 GeV/c^ 
< M{nfi) — mj/^,^,,) < 0.08 GeV/c^. The charmonium 
B DX background can contribute to the muon modes 
if a pion from D meson is misidentified as a muon. Ad- 
ditional veto windows \MKfi — tidI < 0.02 GeV/c^ and 
IMk-jth — mol < 0.02 GeV/c^ suppress this background. 
The invariant mass of the electron pair must exceed 0.14 
GeV/c^ in order to remove background from photon con- 
versions and 7r° — > je~^e~ decays. 

If multiple B candidates survive these selections in an 
event, we select the one with the smallest | Ai?| . The frac- 
tions of multiple B events are about 7%, 12%, and 20% 
for the if+TT", KgTT'^, and A'+7r° modes, respectively, 
according to a Monte Carlo (MC) study. 

We perform an extended unbinned maximum likeli- 
hood fit to Mbc and Mk-k in B — > K*i'^£~ decays, and 
to Mbc in i? — > Ki~^£~ decays. The likelihood function 
is defined as follows: 



hhPK(')hh\ ■ 

where N denotes the number of observed events in the 
candidate region, and TV, (P]), TVt {PI), N,-,x (P^^x). 
and Nj^(,)i^f^ (-P^(.),j^) denote the event yields (the prob- 
ability density functions, PDFs, for the i-th event) for 
signal, combinatorial, B J/'iIj{i/j')X, and B — > K^*^hh 
backgrounds. The signal PDFs consist of a Gaussian 
(Crystal Ball function [8|) in Mbc for the muon (electron) 
modes and a relativistic Breit-Wigner shape in Mktt for 
the K* resonance. The self-cross-feed PDFs, where the 



pion or kaon is misidentified, are modeled by a two- 
dimensional smoothed histogram function and included 
in the signal PDFs as well. The means and widths are de- 
termined from MC and calibrated using 
decays. The combinatorial PDFs are represented by a 
product of an empirical background function introduced 
by ARGUS @ in Mbc, and a threshold function, in which 
the Mktt threshold is fixed at mK+rriTr, plus a relativistic 
Breit-Wigner shape at the K* resonance in Mk-k- The 
PDFs and yields for B J/tp{ip')X decays are deter- 
mined from large MC sample, while the B K^*^hh 
PDFs and normalizations are determined from measured 
data, taking into account the probabilities of the pions 
being misidentified as muons. Yields for signal and com- 
binatorial background, and the combinatorial PDF pa- 
rameters are allowed to float in the fit while the yields 
and parameters for other components are fixed. Fig. [T] 
illustrates the fits for B yields, which are 

230123 and 

166li5 for the K*e+e- and !<£+£- modes, respectively. 




M,, (GeV/c^) 



FIG. 1: Distributions of Mk-k (Mbc) with fit results super- 
imposed for the events in the Mbc (Mktv) signal region. The 
solid curves, solid peak, dashed curves, and dotted curves 
represent the combined fit result, fitted signal, combinatorial 
background, and J/tl;{tp')X background, respectively. 

We divide into 6 bins and extract the signal and 
combinatorial background yields in each bin. The K* 
longitudinal polarization fractions (Pl) and Afb are 
extracted from fits in the signal region to cos 9k* and 
cos 9 Be, respectively, where 9k* is the angle between the 
kaon direction and the direction opposite the B meson 
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in the K* rest frame, and Obi is the angle between the 
and the opposite of the B (B) direction in the 
dilepton rest frame. The signal PDF for the fit to cos 9k* 
is described by 

[ ^Fl cos 0K* + 1(1 - - COS^kO ] X <COS0K') , 

where e{cos9K-') denotes the efficiency obtained from 
MC. For the fit to cosObi, we use 

[ iFUl - cos2 OBt) + f (1 - + cos^ Obi) 

+ApB cosObi ] X e{cos9Bi) 

as the signal PDF, where e(cos0sf) denotes the efficiency 
as a function of cos Obi- The angular efficiency distri- 
butions, background PDFs, and signal and background 
sizes, obtained from cither MC or a Mbc-Mifjr fit, are 
fixed in both angular fits. Fl [Aps] is the only free 
parameter in the fit to cos^if- (cos^s^). Table [J lists 
the measurements of B yields, F^, ApB, and the partial 
branching fractions, obtained by correcting the B yields 
for dependent efficiencies. The differential branching 
fraction, Fl, and Aps as functions of for K*i^£^ and 
Ki+£~ modes are shown in Fig. [2l Fig. [3l and Fig. gl 
respectively. The total branching fractions for the en- 
tire q^ region, obtained by extrapolation from the partial 
branching fractions, as well as the CP asymmetries, for 
the B K*£+£- and B Ki+£- modes are listed in 
Table ini 

We calculate the ratios of branching fractions for the 
electron mode to the muon mode. The lepton flavor ratio 
for B K*£~^£~ [Rk* ) is sensitive to the size of the 
photon pole and is predicted to be 1.33 in the SM, while 
the ratio for B K£~^£~ (Rk) is sensitive to Higgs 
emission and predicted to be larger than 1.0 in the Higgs 
doublet model with large tan/3 The results are 
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FIG. 2: Differential branching fractions for (a) K*£^£^ and 
(b) K£'^£~ modes as a function of q^. The two shaded regions 
are veto windows to reject J/jMtp')X events. The solid curve 
is the theoretical prediction [1^ . 



Rk* = 1.21 ±0.25 ±0.08 , 
Rk = 0.97 ±0.18 ±0.06 . 

Assuming the ratios of branching fractions for the elec- 
tron mode to the muon mode is 1.33 (1.0) in the K*{K) 
mode, the combined branching fractions are measured to 
be 



B{B ^ K*£+£-) = (10.8j:i:J±0.9) X 10"^ , 
B{B^K£+£-) ^ (4.8i°:^±0.3) X 10"^ . 

Isospin asymmetry, shown in TableUand Fig.[5l is defined 



A,= 



_ {tb+Itbo) X B{K^*^°£+£-) - B{K^*'>^£+£-) 
{tb+/tbo) X B{K(*)f^£+£-) + B{Ki*)^£+£-) 




10 12 14 16 18 20 



q^(GeVV) 



FIG. 3: Fit results for Fl as a function of . The solid 
(dashed) curve shows the SM (C7 = —Cf^) prediction. 



where Tb+/tb« (=1.071) is the lifetime ratio of B'^ to 
B^ [llj . A large isospin asymmetry for q^ below the J/V' 



resonance was reported recently [12| . We also measure 
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TABLE I: Fit results in each ol 6 
available [l3|. 



bins and an additional bin from 1 to 6 GeV^ /(? for which recent theory predictions are 



(GeV7c") 



A, 



Fl 



B ^ K* 



6.6 



1.291°: 

0.991J5: 
2.24l°- 
1.051°; 
2.04lij' 
T49^ 



+o.i« 

-0.16 



0.29lt;;i|±0.02 
0-75lo.22±0-05 
0.65l^;2?±0.06 
0-17lo:l5±0.03 
-0.15lg:23±0.07 
0.12lo.l3±0.02 



0.141°^ 
0-47lo: 
0.43l°; 
0-70lo: 
0.66lo 



0-2 
2-5 
5-8.68 
10.09-12.8e 
14.18-16 
>16 



27.4 



25.5lg g 
20.21*:^ 

-10.5 
-9.6 

36.2l^J 

+ 11.0 
-9.9 



35 ±0.12 



-0.67 



±0.10 i.ir 



+0.72 



±0.03 
±0.02 



0.47: 



54.0^ 



84.43 



36 ±0.08 
1^±0.18 
26 ±0-08 
?I±0.16 



-0.47l|^J9±0.04 



0.00 



-8.0 



+0.20 
-0.21 
-+0.30^ 



±0.05 



0-16_o35±0.05 

-,+0.20_ 



0.33l^^^±0.05 



0-67lo,23±0.05 0.26 



32 ±0.03 
26 ±0.07 
25 ±0-14 
2o±0-03 
22 ±0-10 
J«±0.04 



1-6 



29.42 



±0.12 



±0.07 



B -> Kt+r 



-5.4 

22.51^:° 

+7.1 
6.5 
+6.2 
5.5 
+4.9 
4.3 
+8.2 
7.5 



0-8i±o;i6±0-05 

0.58l°;J|±0.04 
0.86l!j;l|±0.05 
0.55lo;i4±0.03 
0-38lo;il±0.02 
0.98lg j^±0.06 



D+0.33_ 

J-a.25- 
-0.49i;jj^±0-04 

-0-19lo.22±0.05 
-0.29lo.29±0.05 
-0.40i;5jg±0.04 
0.1llg:2i±0.05 
"0-4llo.2o±0-04 



0.06l^j5±0.02 
-0.5lt'^|i±0.09 
-0.18t;^ig±0.03 
-0.2lt^lg±0.06 
0-04l:o.26±0.05 
0.02l:°-J,J±0.02 



0-2 
2-5 
5-8.68 
10.09-12.86 
14.18-16 
>16 



27.0; 



34.1 
22.0 
15.6 
40.3 



52.0180 1.36lo;2i±0.08~ 



-0.04lo;^g±0.05 
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TABLE II: Total branching fractions for B K'l+r and 
B K£+r decays. 



Mode 


B (10-') Acp 


K* jj,^ 


11.4l|7±1.0 -0.12l^;^|±0.02 
10.8li;5±0.7 O.OOl!5;i5±0.03 
11.2l^;^±0.8 -0.03lJ5;i3±0.02 


K'+ee 
K'-°ee 
K*ee 


16.4i;^;^±1.8 -0.14l^;^^±0.02 
11.8l2;i±0.9 -0.2llg;;^±0.02 
13.7l2;o±1.2 -0.18lo;i5±0.01 


K'+M 
K*°££ 
K*££ 


12.4l2;o±1.2 -0.13l^;i^±0.01 
9.8ll;i±0.7 -0.08lJ5;i2±0.02 
10.8li;J±0.9 -0.10l!5;i;;±0.01 




5.3lo;7±0.3 -0.05lo;i3±0.03 

4.3lio±0.2 

5.0l°:^±0.3 


K+ee 
K°ee 
Kee 


5.7lo;|±0.3 -0.14l^;J:^±0.03 

2.0li;;^±0.1 

4.8lo;?±0.3 


K+££ 
K°££ 
K££ 


5.3lo;5±0.3 -0.04l^;i^±0.02 

3.3lJ5;?±0.2 

4.8lJ5;^±0.3 



the combined Aj for < 8.68 GeV^/c^ and find 

Ai{B ^ K*l+r) = -0.29l|^:i^±0.03 CT = 1.40, 
Ai{B Ktr) = -0.3ll°;}4±0.05 ct = 1.75 , 



Ai{B ^ K'^*k+r) = -0.30l^:}i±0.04 CT = 2.24, 

where a denotes the significance from null asymmetry 
and is defined as ct = -^Z— 21n {Cq / Cmax) , where Cq is 




10 12 14 16 18 20 



q^(GeVV) 



FIG. 4: Fit results for Afb as a function of q^. The solid 
(dashed) curve shows the SM (C7 — —C^'^') prediction. 



the hkelihood with Aj constrained to be zero and £max 
is the maximum likelihood. Systematic uncertainties are 
considered in the significance calculation. 

Systematic uncertainties in the branching fraction 
measurement for each decay channel are summarized in 
Table Hill They stem dominantly from tracking efficien- 
cies (2.0%-4.4%), MC decay models (0.9%-4.6%), elec- 
tron (3.0%) and muon (2.6%) identification, (4.9%) 
and 7r° (4.0%) reconstruction, and TZ and TZb selec- 
tion (1.2%-3.6%). The signal MC samples are generated 
based on a decay model derived from [l3| , and the mod- 
eling uncertainties are evaluated by comparing different 
MC samples based on different decay models 1^, while 
lepton identification is studied using a, J/tjj ^ £~^i~ data 
control sample. For TZ and TZb selections, we estimate the 
uncertainties from large control samples with the same 
final states, B J/ipK^*'^ with J/ip £+£~. Other 
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TABLE III: Systematic uncertainties in the branching fraction (in %) for each decay channel of electron/muon mode. 



Source 


(T^+TT-)*" 






K+ Ks 


Tracking 


4.4 / 4.4 


3.3 / 3.3 


3.3 / 3.3 


3.0 / 3.0 2.0 / 2.0 


K or 


1.0 / 1.0 


4.9 / 4.9 


1.0 / 1.0 


1.0 / 1.0 4.9 / 4.9 


n or tt" 


1.0 / 1.0 


1.0 / 1.0 


4.0 / 4.0 


- / - - / - 


e/fi 


2.6 / 3.0 


2.6 / 3.0 


2.6 / 3.0 


2.6 / 3.0 2.6 / 3.0 


TZ and TZb 


2.1 / 2.0 


1.9 / 3.6 


3.0 / 3.4 


2.0 / 2.5 2.5 / 1.2 


MC model 


2.5 / 2.9 


2.1 / 2.9 


4.6 / 2.1 


0.9 / 1.0 2.9 / 2.6 


Fitting PDF 


1.6 / 2.1 


3.0 / 3.3 


4.3 / 4.6 


1.4 / 2.0 2.5 / 2.1 


BB pairs 


1.4 / 1.4 


1.4 / 1.4 


1.4 / 1.4 


1.4 / 1.4 1.4 / 1.4 


Rare B 


0.6 / 0.5 


0.9 / 0.9 


0.3 / 0.3 


1.3 / 1.0 1.3 / 0.8 


ccX 


1.1 / 1.4 


1.8 / 2.0 


4.2 / 7.8 


0.3 / 0.4 0.1 / 0.3 


Total 


6.5 / 7.3 


8.0 / 10.1 


9.8 / 12.0 


5.1 / 5.8 7.4 / 7.2 



r- 2 

1.5 
1 

0.5 



-0.5 

-1 - 

-1.5 r 

-2 




2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 
q^(GeVV) 



FIG. 5: Ai as a function of for K*l+r (red) and Ki+T 
(blue) modes. 



TABLE IV: Systematic errors on Fl and Afb measurements. 



Source 


Fl 


Afb 


Signal yield 


0.01-0.06 


0.00-0.06 


Background 


0.01-0.03 


0.01-0.03 


Fl 




0.01-0.13 


Fitting bias 


0.01 


0.02 


Fitting PDF 


0.01 


0.01 



uncertainties such as kaon and pion identification effi- 
ciencies, fitting PDFs, background contamination from 
J /if: decays and charmless B decays, and the number of 
BB pairs are found to be small. The systematic error on 
Rxi") {Ai) is determined by combining the uncertainties 
from lepton (K/n) identification, TZ and TZb selections, 
fitting PDFs and background contamination. The un- 
certainty in Aj from the assumption of equal production 
of i?" and B^ is also considered. Table IIVI shows the 
systematic uncertainties for angular fits. The main un- 
certainties are propagated from the errors on the fixed 
normalizations and Fl, determined from M^jc-Mktt and 



cos9ji* fits, respectively. Fitting bias and fitting PDFs 
are checked using large B — > J/tpK'^*^ and MC samples. 
The total uncertainties range from 0.02-0.06 and 0.03- 
0.15 for Fl and Afb fits, respectively. The systematic 
errors on Ac p are assigned using the measured CP asym- 
metry for sideband data without TZ and TZb selections 
and are found to be around 0.01-0.03. 

In summary, we report the differential branching frac- 
tion and isospin asymmetry as functions of g^, lepton 
flavor ratios, and CP asymmetries in both B K*t^l~ 
and B K(,^£~ decays. K* longitudinal polarization 
and forward-backward asymmetry as functions of in 
B K*£'^£~ are also measured from an angular analy- 
sis. The differential branching fraction, lepton flavor ra- 
tios, and K* polarization are in good agreement with the 
SM predictions. No significant CP asymmetry or isospin 
asymmetry is found. The Apb{<i'^) spectrum, although 
consistent with previous measurements [lij . tends to be 
shifted toward the positive side from the SM expecta- 
tion, especially at large values. A much larger data 
set, such as will be available from the proposed super 
B factory [31 and LHCb 17|, is needed to make more 



precise comparisons with the SM and other theoretical 
predictions. 

We thank the KEKB group for the excellent opera- 
tion of the accelerator, the KEK cryogenics group for the 
efficient operation of the solenoid, and the KEK com- 
puter group and the National Institute of Informatics 
for valuable computing and SINET3 network support. 
We acknowledge support from the Ministry of Educa- 
tion, Culture, Sports, Science, and Technology of Japan 
and the Japan Society for the Promotion of Science; the 
Australian Research Council and the Australian Depart- 
ment of Education, Science and Training; the National 
Natural Science Foundation of China under Contracts 
No. 10575109 and 10775142; the Department of Science 
and Technology of India; the BK21 program of the Min- 
istry of Education of Korea, the CHEP SRC program 



8 



and Basic Research program (Grant No. ROl-2005-000- 
10089-0) of the Korea Science and Engineering Founda- 
tion, and the Pure Basic Research Group program of the 
Korea Research Foundation; the Pohsh State Committee 
for Scientific Research; the Ministry of Education and 
Science of the Russian Federation and the Russian Fed- 
eral Agency for Atomic Energy; the Slovenian Research 
Agency; the Swiss National Science Foundation; the Na- 
tional Science Council and the Ministry of Education of 
Taiwan; and the U.S. Department of Energy. 



[1] A. All, P. Ball, L. T. Handoko and G. Hiller, Phys. Rev. 
D 61, 074024 (2000); Y. Wang and D. Atwood, Phys. 
Rev. D 68, 094016 (2003). 

[2] S. Davidson, D. C. Bailey and B. A. Campbell, Z. Phys. 
C 61, 613 (1994); A. Ah, G. F. Giudice and T. Mannel, 
Z. Phys. C 67, 417 (1995); G. Burdman, Phys. Rev. D 
52, 6400 (1995); J. L. Hewett and J. D. Wells, Phys. 
Rev. D 55, 5549 (1997); A. Ah, E. Lunghi, C. Greub and 
G. Hiller, Phys. Rev. D 66, 034002 (2002); T. M. Aliev, 
A. Ozpineci and M. Savci, Eur. Phys. J. C 29, 265 (2003) 
and references therein. W. S. Hon, A. Hovhannisyan and 
N. Mahajan, Phys. Rev. D 77, 014016 (2008). 

[3] S. Kurokawa and E. Kikutani, Nucl. Instrum. Methods 
Phys. Res., Sect. A 499, 1 (2003) and other papers in- 
cluded in this Volume. 

[4] A. Abashian et al. (Belle Collaboration), Nucl. Instrum. 
Methods Phys. Res., Sect. A 479, 117 (2002). 

[5] K. Hanagaki et al, Nucl. Instrum. Methods Phys. Res., 
Sect. A 485, 490 (2002); A. Abashian et al, Nucl. In- 
strum. Methods Phys. Res., Sect. A 491, 69 (2002). 

[6] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 
(1978). The modified moments used in this paper are 
described in S. H. Lee et al. (Belle Collaboration), Phys. 



Rev. Lett. 91, 261801 (2003). 
[7] H. Kakuno et al, Nucl. Instrum. Methods Phys. Res., 

Sect. A 533, 516 (2004). 
[8] T. Skwarnicki (Crystal Ball Collaboration), Ph.D. thesis, 

Cracow Institute of Nuclear Physics, DESY F31-86-02 

(1986). 

[9] H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. 

B 185, 218 (1987). 
[10] C. Bobeth, G. Hiller and G. Piranishvili, 

larXiv:0709.4T74l M. Beneke, Th. Feldmann and 

D. Seidel, Eur. Phys. J. C 41, 173 (2005). 
[11] C. Amsler et al. (Particle Data Group), Phys. Lett. B 

667, 1 (2008). 

[12] B. Aubert et al. (BaBar Collaboration), arXiv:hep- 
ex/0807.4119 (2008). 

[13] A. Ah, E. Lunghi, C. Greub and G. Hiller, Phys. Rev. D 
66, 034002 (2002). 

[14] D. Melikhov, N. Nikitin and S. Simula, Phys. Lett. B 
410, 290 (1997); P. Colangelo, F. DeFazio, P. Santorelh 
and E. Scrimieri, Phys. Rev. D 53, 3672 (1996). 

[15] K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 
88, 021801 (2001); A. Ishikawa et al. (Belle Collabora- 
tion), Phys. Rev. Lett. 91, 261601 (2003); M. Iwasaki 
et al. (Belle Collaboration), Phys. Rev. D 72, 092005 
(2005); A. Ishikawa et al. (Belle Collaboration), Phys. 
Rev. Lett. 96, 251801 (2006); B. Aubert et al. (BaBar 
Collaboration), Phys. Rev. Lett. 93, 081802 (2004); 
B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 
73, 092001 (2006); B. Aubert et al. (BaBar Collabora- 
tion), arXiv:hep-ex/0804.4412vl (2008). 

[16] A. G. A keroyd et al. (The Sup erKEKB Physics Working 
Group) , |arXiv:hep-ex/0406071| 

[17] J. Dickens (LHCb Collaboration), in Proceedings of the 
CKM 2006 Workshop, Nagoy a, Japan, 200 6 (unpub- 
lished, http://agenda.hepl.phys.nagoya-u.ac.jp/g etFile] 
py / access?contribId= 18&sessionld=9&:resld=0& 
materialId=slides&confId=6 . ) 



